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Abstract 



The ALICE experiment has measured the inclusive production in Pb-Pb collisions at -^Inn 



2.76 TeV down to pt = in the rapidity range 2.5 < y < 4. A suppression of the inclusive 
yield in Pb-Pb is observed with respect to the one measured in pp collisions scaled by the number 
of binary nucleon-nucleon collisions. The nuclear modification factor, integrated over the 0%-80% 
most central collisions, is 0.545 ± 0.032(stat.) ± 0.084(syst.) and does not exhibit a significant de- 
pendence on the collision centrality. These features appear significantly different from lower energy 
measurements. Models including J/xj/ production from charm quarks in a deconfined partonic phase 
can describe our data. 



*See Appendix lAlfor the list of collaboration members 
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Ultra-relativistic collisions of heavy nuclei aim at producing nuclear matter at high temperature and 
pressure. Under such conditions Quantum Chromodynamics predicts the existence of a deconfined state 
of partonic matter, the Quark-Gluon Plasma (QGP). Among the possible probes of the QGP, heavy quarks 
are of particular interest since they are expected to be produced in the primary partonic scatterings and 
to coexist with the surrounding medium. Therefore, the measurement of quarkonium states and hadrons 
with open heavy flavor is expected to provide essential information on the properties of the strongly- 
interacting system formed in the early stages of heavy-ion collisions HI. In particular, according to the 
color-screening model El, measuring the in-medium dissociation probability of the different quarkonium 
states is expected to provide an estimate of the initial temperature of the system. In the past two decades, 
J/\j/ production in heavy-ion collisions was intensively studied at the Super Proton Synchrotron (SPS) 
and at the Relativistic Heavy Ion CoUider (RHIC), from ^/s^ « 20 to 200 GeV. At the SPS, a strong J/}]/ 
suppression was found in the most central Pb-Pb collisions [3|. The observed suppression is larger than 
the one expected from Cold Nuclear Matter (CNM) effects, which include nuclear absorption and (anti-) 
shadowing. The dissociation of excited cc states like Xc and i//^(2S), which in pp collisions constitute 
about 40% of the inclusive yield, is one possible interpretation of the observed suppression. A 
J/\j/ suppression was also observed at RHIC, in central Au-Au collisions HQ, at a level similar to 
the one observed at the SPS when measured at mid-rapidity although it is larger at forward rapidity. 
Several models |7l [8] |9l attempt to reproduce the RHIC data by adding to the direct J/i/a production 
a regeneration component from deconfined charm quarks in the medium, which counteracts the J/i/a 
dissociation in a QGP. A quantitative description of these final-state effects is however difficult at the 
present time because of the lack of precision in the CNM effects and in the open charm cross section 
determination. The measurement of charmonium production is especially promising at the Large Hadron 
Collider (LHC) where the high energy density of the medium and the large number of cc pairs produced 
in central Pb-Pb collisions should help to disentangle between the different suppression and regeneration 
scenarios. At the LHC, a suppression of inclusive J/i// with high transverse momentum was observed in 
central Pb-Pb collisions at t/5nn = 2.76 TeV with respect to peripheral collisions or pp collisions at the 
same energy by ATLAS UHl and CMS HD, respectively. 

In this Letter, we report ALICE results on inclusive J/\j/ production in Pb-Pb collisions at ^/5nn = 2.76 
TeV at forward rapidity, measured via the H^jJ.^ decay channel. Our measurement encloses the low 
transverse momentum region that is not accessible to other LHC experiments and thus complements 
their observations. The J/i/a corrected yield in Pb-Pb collisions is combined with the one measured in pp 
collisions at the same center-of-mass energy |[T2l to form the J/i/a nuclear modification factor Raa- The 
results are presented as a function of collision centrality and rapidity (y), and in intervals of transverse 
momentum (pi). 

The ALICE detector is described in |[T3]| . At forward rapidity (2.5 <y <4) the production of quarkonium 
states is measured in the muon spectrometer Q down to pt = 0. The spectrometer consists of a ten 
interaction length thick absorber filtering the muons in front of five tracking stations comprising two 
planes of cathode pad chambers each, with the third station inside a dipole magnet with a 3 Tm field 
integral. The tracking apparatus is completed by a triggering system made of four planes of resistive plate 
chambers downstream of a 1.2 m thick iron wall, which absorbs secondary hadrons escaping from the 
front absorber and low momentum muons coming mainly from n and K decays. In addition, the silicon 
pixel detector (SPD) and scintillator arrays (VZERO) were used in this analysis. The VZERO counters, 
two arrays of 32 scintillator tiles each, cover 2.8 < T] <5.1 (VZERO-A) and -3.7 < T] < -1.7 (VZERO- 
C). The SPD consists of two cylindrical layers covering |tj| < 2.0 and |t]| < 1.4 for the inner and outer 
layers, respectively. All these detectors have full azimuthal coverage. The minimum bias (MB) trigger 
requirement used for this analysis consists of a logical AND of the three following conditions: (i) a signal 
in two readout chips in the outer layer of the SPD, (ii) a signal in VZERO-A, (iii) a signal in VZERO- 



In the ALICE reference frame, the muon spectrometer covers a negative rj range and consequently a negative y range. We 
have chosen to present our results with a positive y notation. 
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C, providing a high triggering efficiency for hadronic interactions. The beam induced background was 
further reduced by timing cuts on the signals from the VZERO and from the zero degree calorimeters 
(ZDC). The contribution from electromagnetic processes was removed with a cut on the energy deposited 
in the neutron ZDCs. The centrality determination is based on a fit of the VZERO amplitude distribution 
as described in [14|. A cut corresponding to the most central 80% of the nuclear cross section was 
applied; for these events the MB trigger is fully efficient. A data sample of 17.7 x 10^ Pb-Pb collisions 
collected in 2010 satisfying all the above conditions is used in the following analysis. It corresponds 
to an integrated luminosity = 2.88 ±0.10 /ib^'. This data sample was further divided into five 
centrality classes from 0%-10% (central collisions) to 50%-80% (peripheral collisions). 

J/Y candidates are formed by combining pairs of opposite-sign (OS) tracks reconstructed in the geomet- 
rical acceptance of the muon spectrometer. To reduce the combinatorial background, the reconstructed 
tracks in the muon tracking chambers are required to match a track segment in the muon trigger sys- 
tem. The J/i/A raw yield was extracted by using two different methods. The OS dimuon invariant mass 
distribution was fitted with a Crystal Ball (CB) function, i.e. a Gaussian with a polynomial tail ifTSll . 
to reproduce the J/y/ line shape, and a sum of two exponentials to describe the underlying continuum. 
At high transverse momenta (pt > 3 GeV/c), the sum of two exponentials does not describe correctly 
the underlying continuum; it was replaced by a third order polynomial. Alternatively, the combinatorial 
background was subtracted using an event-mixing technique. The resulting mass distribution was fitted 
with a CB function and an exponential or a first order polynomial to describe the remaining background. 
The V'^(2S) was not included in the signal line shape since its contribution is negligible. The width of 
the J/i/A mass peak depends on the resolution of the spectrometer which can be affected by the detector 
occupancy that increases with centrality. This effect, evaluated by embedding simulated J/i/A — )• 
decays into real events, was found to be less than 2%. This conclusion was confirmed by a direct mea- 
surement of the tracking chamber resolution versus centrality using reconstructed tracks. Therefore, the 
same CB line shape can be used for all centrality classes. The parameters of the CB tails were fixed to 
the values obtained either in simulations or in pp collisions where the signal to background ratio is much 
higher. For each of these choices, the mean and width of the CB Gaussian part were fixed to the value 
obtained by fitting the mass distribution in the centrality range 0%-80%. In addition, a variation of the 
width of lb 1 standard deviation was applied to account for uncertainties (varying the mean has turned 
out to have a negligible effect in comparison). The raw J/i/a yield in each centrality class was determined 
as the average of the results obtained with the two fitting approaches and the various CB parametriza- 
tions, while the corresponding systematic uncertainties were defined as the RMS of these results. It was 
also checked that every individual result differs from the mean value by less than three RMS. The raw 
J/i//^ yield in our Pb-Pb sample is 2350± 139(stat.) it 189(syst.). The invariant mass resolution is around 
75 MeV/c^, in very good agreement with the embedded J/y simulations. The signal to background ratio 
integrated over zt 3 a of the mass resolution varies from 0.1 for central collisions to 1.5 for peripheral 
colUsions. 

The measured number of J/y (^j/^,) was normalized to the number of events in the corresponding cen- 
trality class (A'^vents) ^^'^ further corrected for the branching ratio (BR) of the dimuon decay channel, the 
acceptance A and the efficiency e' of the detector. The inclusive J/y yield in each centrality class / for 
our measured pi and y ranges {Api, Ay) is then given by: 

N', 

YLJApuAy) = :. (1) 

^ ^J/ + events^ 

The product Ae was determined from Monte Caiio simulations. The generated J/y Pi and y distributions 
were extrapolated from existing measurements [16], including shadowing effects from EKS98 calcula- 
tions [17|. As J/y mesons produced from charm quarks in the medium are expected to be unpolarized, 
and the measured J/y polarization in pp collisions at = 7 TeV is compatible with zero ifTSl . we 
presume J/y production is unpolarized. For the tracking chambers, the time-dependent status of each 
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Table 1: The average number of binary nucleon-nucleon collisions (A^coii). of participating nucleons (A^part) without 
and with (A^coii) weighting, the mid-rapidity charged-particle density d^V^/dTjjjj^o with (A^cou) weighting and the 
average value of the nuclear overlap function (Taa) for the centrality classes expressed in percentages of the 
nuclear cross section lfT4l . 



Centrality 


(A^coll) 


(A^part) 


(A^p^art) 


di7lD=o 


{Taa) (mb-i) 


0%-10% 


1503±170 


356±4 


361 ±4 


1463±60 


23.5±1.0 


10%-20% 


923±100 


260±4 


264±4 


979±37 


14.4±0.6 


20%-30% 


559±56 


186±4 


189±4 


658±23 


8.74±0.37 


30%-50% 


248 ±24 


107±3 


117±3 


369±13 


3.87±0.18 


50%-80% 


46±4 


32±2 


47 ±2 


110±5 


0.72±0.05 


0%-80% 


450±48 


139±3 


264±4 


n.a. 


7.03±0.27 



electronic channel during the data taking period was taken into account as well as the residual misalign- 
95 ment of the detection elements. The efficiencies of the muon trigger chambers were determined from 
data and were then applied in the simulations [ 19 1. Finally, the dependence of the efficiency on the de- 
tector occupancy was included using the embedding technique. For J/i/a mesons emitted at 2.5 <y < 4 
and pi>0,a. run-averaged value of Ae = 0.176 with a 8% relative systematic uncertainty was obtained. 
A 8% ±2%(syst.) relative decrease of the efficiency was observed when going from peripheral to central 
100 collisions. 

The J/i/A yield measured in Pb-Pb collisions in centrality class / is combined with the inclusive J/i/a cross 
section measured in pp collisions at the same energy to form the nuclear modification factor Raa defined 
as: 



The inclusive J/i/a cross section in pp collisions (JjJ^{Apt,Ay) was measured using the same apparatus and 
105 analysis technique within the corresponding pt and y range [ 12|. The reference value a^J^ used for the 

calculation of /?aa integrated over pt andy is 3.34±0.13(stat.) ±0.24(syst.) ±0.14(lumi.)+'j'o7(pol.) /ib. 
The centrality intervals used in this analysis and the corresponding average number of binary collisions 
(Ncou), average number of participating nucleons (A'part) and average value of the nuclear overlap func- 
tion (7aa) derived from a Glauber model calculation [14] are summarized in Table [T] Since our most 

no peripheral bin is rather large, the variables (A'^rt) and the charged-particle density measured at mid- 
rapidity dN^^^/dri\jj=o were weighted by the number of binary collisions (A'coii)- Indeed in absence of 
nuclear matter effects, the J/y production cross section in A- A is expected to scale with (A'coii)- The 
weighted values are given in Table[T]and are used for the ALICE data points in the following figures. All 
systematic uncertainties entering the Raa calculation are listed in Table |2l In the figures below, the point 

115 to point uncorrected systematic uncertainties are represented as boxes at the position of the data points 
while the statistical ones are indicated by vertical bars. Correlated systematic uncertainties are quoted 
directly on the figures. 



The inclusive J/i/A Raa measured by ALICE at y^^NN = 2.76 TeV in the range 2.5 < y <4 and pt>0 
is shown in Fig. [T] as a function of dA'^ch/dr/ |rj=o (left) and A'part (right). The charged-particle density 
120 closely relates to the energy density of the created medium whereas the number of participant reflects the 
collision geometry. The centrality integrated J/xj/ Raa is /?aa"^°'* = 0.545 ± 0.032(stat.) ± 0.084(syst.) , 
indicating a clear J/y suppression. The contribution from beauty hadron feed-down to the inclusive J/y 
yield in our y and pt domain was measured by the LHCb collaboration to be about 10% in pp colli- 
sions at ^/s = 7 TeV 1211. Therefore, the difference between the prompt J/i/a /?aa and our inclusive 
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Table 2: Summary of the systematic uncertainties entering the /?aa calculation. The type I (II) stands for correlated 
(uncorrelated) uncertainties. The centrality dependence for the type II is given as a range. 



source value type 

signal extraction 5%-12% II 

input MC parametrization 5% I 

tracking efficiency 5%andO%-l% I and II 

trigger efficiency 4% and 0%-2% I and II 

matching efficiency 2% I 

Taa 4%-8% n 

ofL at = 2.76 TeV 9% I 



< i-4r 

of : 

1.2- 



■ ALICE (Pb-Pb \js^ = 2.76 TeV), 2.5<y<4 

• PHENIX (Au-Au ^ = 200 GeV), 1 .2<|y|<2.2 

O PHENIX (Au-Au \ls;;^ = 200 GeV), |y|<0.35 




global sys.= ± 1 2.5% 
global sys.= ± 9.2% 
global sys.= ± 1 2% 



200 400 600 800 1000 1200 1400 



m ALICE (Pb-Pb ^ = 2.76 TeV), 2.5<y<4 global sys.= ± 12.5% 
• PHENIX (Au-Au ^ = 200 GeV), 1 .2<|y|<2.2 global sys.= ± 9.2% 
o PHENIX (Au-Au ^= 200 GeV), |y|<0.35 global sys.= + 12%, 
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Fig. 1: (Color online) Inclusive J/i/a R^a as a function of the mid-rapidity charged-particle density (left) and the 
number of participating nucleons (right) measured in Pb-Pb collisions at ^/Imn = 2.76 TeV compared to PHENIX 
results in Au-Au collisions at ^/s^ = 200 GeV at mid-rapidity and forward rapidity 1.4. ,5. ,20J . 
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Fig. 2: (Color online) Centrality integrated (0%-80%) inclusive J/v//?aa rneasured in Pb-Pb collisions at -^^nn = 
2.76 TeV as a function of rapidity for two pt ranges. The open boxes contain the total systematic uncertainties 
except the one from the integrated luminosity in the pp reference of 4.3% (6%) for the ALICE (CMS |1JJ) data. 
The two models ||22l|23] predict the 7?aa due only to shadowing effects for nDSg (shaded areas) and EPS09 (Unes) 
nPDF respectively. 



125 measurement is expected not to exceed 11% if Ncoi\ scaling of beauty production is assumed and shad- 
owing effects are neglected. The comparison with the PHENIX measurements H at y^^NN = 200 GeV 
at forward rapidity 1.2 < j^^l < 2.2 ||5l|23 shows that our inclusive J/y Raa is almost a factor of three 
larger for dA^ch/dT] |,j=o ^ 600 (A'part > 180). In addition, our results do not exhibit a significant centrality 
dependence. 

130 The rapidity dependence of the J/i// Raa is presented in Fig. |2] for two pi domains, pt> and pi > 
3 GeV/c. The J/ i/A reference cross sections in pp collisions and the Raa total systematic uncertainties, 
indicated as open boxes in the figure, were evaluated in the same kinematic range. Our results are 
shown together with a measurement from CMS 1 1 1 1 of the inclusive J/i/a Raa in the rapidity range 
1.6 < \y\ < 2.4 with pt>3 GeV/c. No significant rapidity dependence can be seen in the J/\j/ Raa 

135 for Pt > 0. For pt>3 GeV /c, a decrease of Raa is observed with increasing rapidity reaching a value 
of 0.289 ±0.061 (stat.) ±0.078(syst.) for 3.25 <y <4. At LHC energies, J/y nuclear absorption is 
likely to be negligible and the modification of the gluon distribution function is dominated by shadowing 
effects Il24ll . An estimate of shadowing effects is shown in Fig. [2] within the Color Singlet Model at 
Leading Order f22l and the Color Evaporation Model at Next to Leading Order [23]. The shadowing is 

140 respectively calculated with the nDSg and the EPS09 parametrizations of the nuclear Parton Distribution 
Function (nPDF). For nDSg (EPS09) the upper and lower limits correspond to the uncertainty in the 
factorization scale (uncertainty of the nPDF). The effect of shadowing shows no dependence with rapidity 
and its overall amount is reduced by the addition of a transverse momentum cut. At most, shadowing 
effects are expected to lower the Raa from 1 to 0.7. Recent Color Glass Condensate (CGC) calculations 

145 for LHC energies may indicate a larger initial state suppression {Raa ^ 0.5) 1,251 . However, any J/i/A 
suppression due to initial state effects, CGC or shadowing, will be stronger at lower pt contrary to the 
data behavior. 

In Fig. [3l our measurement is compared with theoretical models that include a 3/\j/ regeneration compo- 
nent from deconfined charm quarks in the medium. The Statistical Hadronization Model ||6l|26l assumes 
150 deconfinement and a thermal equilibration of the bulk of the cc pairs. Then charmonium production 
occurs at phase boundary by statistical hadronization of charm quarks. The prediction is given for 

^ The PHENIX mid-rapidity J/y Rw was measured in centrality classes wider than the ones in which the mid-rapidity 
charged-particle density is given 1201 . Therefore a linear interpolation was done to extract the mid-rapidity charged-particle 
density in the three most peripheral classes. 
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Fig. 3: (Color online) Inclusive J/y/ 7?aa measured in Pb-Pb collisions at y'SNN = 2.76 TeV compared to the 
predictions by Statistical Hadronization Model Il26l . Transport Model I fZTl and II f2Fl, see text for details. 



two values of dOcc/dy in absence of a measurement for Pb-Pb collisions. The two transport model 
results ll27l |28]| presented in the same figure differ mostly in the rate equation controlling the J/y/ dis- 
sociation and regeneration. Both are shown as a band which connects the results obtained with (lower 
155 limit) and without (higher limit) shadowing. The width of the band can be interpreted as the uncertainty 
of the prediction. In both transport models, the amount of regenerated J/y in the most central collisions 
contributes to about 50% of the measured yield, the rest being from initial production. 

In summary, we have presented the first measurement of inclusive J/y Raa down to pi = at forward 
rapidity in Pb-Pb collisions at ^/s^m = 2.76 TeV. The J/y Raa is larger than the one measured at the 

160 SPS and at RHIC for most central collisions and does not exhibit a significant centrality dependence. 
Statistical hadronization and transport models which feature a full or partial J/y production from charm 
quarks in the QGP phase can describe the data. Towards a definitive conclusion about the role of J/y 
production from deconfined charm quarks in a partonic phase, the amount of shadowing needs to be 
measured precisely in pPb collisions. In this context, the measurement of open charm and J/y elliptic 

165 flow will also help to determine the degree of thermalization for charm quarks. 
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